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Abstract 

We present a measurement of the decay amplitudes and triple-product correlations in B <pK* 
decays based on 140 fb~^ of data recorded at the T(4S') resonance with the Belle detector at 
the KEKB e+e" storage ring. The decay amplitudes for the different ^K* helicity states are 
measured from the angular distributions of final state particles in the transversity basis. The 
longitudinal and transverse complex amplitude moduli and angles are \ Aq\'^ = 0.51 it 0.06 it 0.04, 
|^_Lp = 0.24±0.06±0.03, arg(^ii) = -2.21 ±0.22±0.05 rad, and arg(^_L) = 0.72±0.21±0.06 rad. 
The T- violating asymmetries through triple-product correlations are measured to be consistent with 
zero. 

PACS numbers: 13.25.Hw, 11.30.Er 
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The decay B (j)K* involves the FCNC 6 s transition that is forbidden to first order 
in the Standard Model (SM) but can proceed by second order penguin and box diagrams. 
This process therefore provides information on the Cabibbo-Kobayashi-Maskawa matrix el- 
ement Vts (H and is sensitive to physics beyond the SM. The decay B — > (pK* is a mixture 
of CP-even and CP-odd states, and polarization measurements allow us to project out the 
different CP states statistically. Measurements of the T-odd triple-product correlations pro- 
vide information about T-violating asymmetries |2| as well as observables that are sensitive 
to New Physics (NP) j^. Our previous measurement 0] and a recent report by BaBar 
both suggest that the fraction of longitudinal polarization component differs from the naive 
SM prediction. This provides a strong motivation to update our study with a larger dataset. 

In this paper, we present a measurement of the decay amplitudes in P+ —>■ (pK*^ and 
P° — > (f)K*^ decays by a full three-dimensional angular analysis. The charge conjugate 
modes are included everywhere unless otherwise specified. We also report the measurements 
of triple-product correlations and related T-violating asymmetries, as well as the observables 
that are sensitive to NP. 

This analysis is based on a data set with an integrated luminosity of 140 fb~^ taken at 
the T(4S') resonance recorded by the Belle detector [6J at the KEKB e~^e~ collider 0. This 
luminosity corresponds to 152 x 10^ produced BB pairs. The beam energies are 8 GeV for 
e~ and 3.5 GeV for e"*". 

The Belle detector is a general purpose magnetic spectrometer equipped with a 1.5 T su- 
perconducting solenoid magnet. Charged tracks are reconstructed in a Central Drift Cham- 
ber (CDC) and a Silicon Vertex Detector (SVD). Photons and electrons are identified using 
a CsI(Tl) Electromagnetic Calorimeter (ECL) located inside the magnet coil. Charged par- 
ticles are identified using specific ionization (dE/dx) measurements in the CDC as well as 
information from Aerogel Cherenkov Counters (ACC) and Time of Flight Counters (TOP). 
A kaon likelihood ratio, P{K/n) = LkI{C.k + C.-k)-, has values between (likely to be a 
pion) and 1 (likely to be a kaon), where C.k{-k) is derived from the dE/dx, ACC and TOP 
measurements. 

Candidate — >• K^K~ decays are found by selecting pairs of oppositely charged tracks 
that are not pion-like {P{K/7r) > 0.1). The vertex of the candidate charged tracks is 
required to be consistent with the interaction point (IP) to remove poorly measured tracks. 
In addition, candidates are required to have a K^K~ invariant mass that is less than 10 
MeV/c^ from the nominal meson mass. 

Pairs of oppositely charged tracks are used to reconstruct — >• tt^tt^ decays. The 7r"'"7r~ 
vertex is required to be displaced from the IP by a minimum transverse distance of 0.22 cm for 
high momentum (> 1.5 GeV/c) candidates and 0.08 cm for those with momentum less than 
1.5 GeV/c. The direction of the pion pair momentum must agree with the direction defined 
by the IP and the vertex displacement within 0.03 rad for high-momentum candidates, and 
within 0.1 rad for the remaining candidates. 

Charged tracks with P{K/7t) > 0.4 (< 0.9) are considered to be kaons (pions). For 
7r° — s> 77, a minimum photon energy of 50 MeV is required and the 77 invariant mass must 
be less than 16 MeV/c^ from the nominal vr*^ mass. K* candidates are reconstructed in three 
decay modes: K*^ — > K^tt~, K*~^ — > K~^Tr^ and K*~^ —>■ Kgii'^. The invariant mass of the 
K* candidate is required to be less than 70 MeV/c^ from the nominal K* mass. 

A B meson is reconstructed from a pair of (p and K* candidates and identified by 
the energy difference AE = E]^ — E^^^^, and the beam constrained mass Mbc = 

\J{E^^^^)^ — {p^B^^y. E^^^^ is the beam energy in the center-of-mass system (cms) of the 
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T(4S') resonance, and i?^™** and p^** are the cms energy and momentum of the reconstructed 
B candidate. The 5-meson signal window is defined as 5.27 GeV/c^ < Mbc < 5.29 GeV/c^ 
and \AE\ < 45 MeV. The signal window is enlarged to -100 MeV < AE < 80 MeV for 
— s> (f)K*^{K*^ K^TT^) because of the impact of shower leakage on AE resolution. An 
additional requirement cos6'x* < 0.8 is applied to reduce low momentum 7r° background, 
where 6k* is the angle between the K* direction and its daughter kaon defined in the K* 
rest frame. In the signal window, about 1% of the events have multiple candidates. We 
choose the best candidate according to the goodness of fit of the B meson's vertex. 

The dominant background is continuum e+e~ —>■ qq production (g = u,d,c,s). Several 
variables are used to exploit the differences between the event shapes for continuum qq 
production (jet-like) and for B decay (spherical) in the cms frame of the T(4S') These 
variables are combined into a single likelihood ratio Rs = Cs/{Cs + Cqq)-, where Cg (Cqq) 
denotes the signal (continuum) likelihood. The requirements on the likelihood ratio are 
determined by a figure of merit study that includes a dependence on the S-tagging quality 

i. 

Backgrounds from other B decay modes such as -B — ^ KKK*, B — > /o(980)K*(/o 
K~^K~), B (pKn, B — >■ KKKn, and feed-across between (pK* and (pK decay channels 
are studied. The contributions from B KKK* and B /o(980)ir*(/o ^ K+K-) are 
estimated from the K^K~ invariant mass distribution. The K^K~~ mass distribution for 
B — > KKK* is determined by Monte Carlo (MC) simulation assuming three-body phase 
space decay. The shape for /o(980) is obtained from MC, where a S'-wave Breit-Wigner with 
a 40 MeV/c^ intrinsic width is assumed These contributions are estimated separately 
by fits to the events outside of the mass region. The contribution from B KKK* is 
estimated to be 5-10% of the signal yield depending on K* decay modes. The B foK* 
contribution is estimated to be 4 — 11%. Contamination from four-body B — > KKKtt decays 
is checked by performing fits to the non-resonant region of K^K~ and Kti mass. It is found 
to be very small and is neglected. To remove the contamination of (j)K decays, they are 
explicitly reconstructed and removed from (pK* candidates. 

The signal yields (A^s) are extracted by extended unbinned maximum-likelihood fits per- 
formed simultaneously in AE and M^^. The signal probability density functions (PDFs) are 
products of Gaussians in AE and Mbc- The means and widths are verified using B ^ J / ipK* 
decays. Additional bifurcated Gaussians (Gaussians with different widths on either side of 
the mean) are used to model the tails in the AE distribution. 

The PDF shape for the continuum events are parameterized by an ARGUS function in 
Mbc and a linear function in AE. The parameters of the functions are determined by a fit 
to the events in the signal sideband. The sideband region is defined by Mbc < 5.265 GeV/c^ 
and \AE\ > 0.08 GeV {AE > 0.10 GeV and AE < -0.12 GeV) for the modes without 
(with) a 7r° in the final state, respectively. The signal and background yields are allowed 
to fioat in the fit while other PDF parameters are fixed. The measured signal yields and 
purities are summarized in Table EJ The distributions of AE and Mbc are shown in Fig. ^ 

The decay angles of a S-meson, to the two vector mesons and K* defined in the 
transversity basis , are shown in Fig. |21 The x-y plane is defined to be the decay plane 
of K* and the x axis is in the direction of the (/)-meson. The y axis is perpendicular to the x 
axis and is on the same side as the kaon from the K* decay. The z axis is perpendicular to 
the x-y plane according to the right-hand rule, 6tr (0tr) is the polar (azimuthal) angle with 
respect to the z-axis of the K^ from (p decay in the rest frame, and 6k* is defined above. 
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FIG. 1: Distributions of AE (Mbc) with fit results for the events in the M^c (AE) signal window. 
The continuum background component is shown by the dashed curves. 



TABLE I: Number of events observed in the signal box {N^y), signal yields {Ns) obtained by fits 
after KKK* background subtraction, and signal purities. 



Mode 








Ns 


Purity 


B^^ 






167 






B+ ^ 
B+ ^ 






40 
52 


20.4+^-^ 





The distribution of decays in the three angles 1J|, 6k*-, ^tr, and 0tr is 

rf3r(0to,COS^tr,COS^X*) 9 2o 2n ■ 2n 



d^xxd cos COS Ok* 327r 
+ l^ii p2 sin^ sin^ 0ti. sin^ Qj^* 

+ |v4o|^4 Sin^ 6'tr COS^ 0tr cos^ Qk* 

+ V2Re(Aj|Ao) sin^ 6'tr sin 20tr sin IQk* 
-r/y2Im(A*Ax)sin20tr cos 0tr sin IQk* 

- 2r]\m.{A\A^_) sin 26^^ sin 0tr sin^ 6'x*] , (1) 

where Aq, Ay, and are the complex amplitudes of the three helicity states in the transver- 
sity basis with the normalization condition \Aq\^ + \A\\^ + |v4_|_p = 1, and r] = +1 (— 1) 
corresponds to B (B) mesons. Ao denotes the longitudinal polarization of the K'^K' 
system and A_[_ {A\\) is the transverse polarization along the 2;-axis (y-axis). The value of 
Ui P (l^oP + l^iiP = 1 - |^±n is the CP-odd (CP-even) fraction in the decay B (pK* 
|l4| . The imaginary phases of the amplitudes are sensitive to final state interactions (FSI). 
The presence of FSI results in phases that differ from either or ±7r. 

The complex amplitudes are determined by performing an unbinned maximum likelihood 
fit with the B (pK* candidates in the signal window. The combined likelihood is given 
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FIG. 2: The definition of decay angles in B ^ (pK* decays. 

by 

N 

i 

+ fgq ■ Rqq{OK*, ^^tr, 0tr) 

+ JkKK* ■ RkKK'^Ok* , dtr, 0tr)] , (2) 

where T is the angular distribution function (ADF) given by Eq. HI and Rqq and Rrkk* 
are the ADFs for continuum and B — > KKK* background, respectively. The value of t] is 
determined from the charge of the kaon in the K* decay, Rqq is determined from sideband 
data and Rrkk* is assumed to be fiat. The detection efficiency function (e) is determined 
by MC. The fractions of (pK* (f^f^K*), ifqq) and B — > KKK* {/kkk*) are parameterized 
as a function of AE and Mbc- The value of arg(Ao) is set to zero and l^yP is calculated 
from the normalization constraint in the fit. The four parameters (l^oP, I^^P; arg(y4||), 
and arg(A^)) are determined from the fit. The results from B^ (pK*^ and B^ — > (pK*^ 
may be combined if both decays are dominated by 6 — s penguin transitions and the 
h ^ u annihilation contribution to B^ cpK*^ can be neglected. Figure El shows the angle 
distributions with projections of the fit superimposed. The amplitudes obtained from the 
fit are shown in Table |nj 

The systematic uncertainties include the slow pion detection efficiency (6 — 7%), the 
background from higher K* states (2 — 8%), and the B — >■ foK* background (2 — 3%). The 
systematic uncertainty in the angular resolution is estimated by MC simulation and found to 
be less than 1%. Uncertainties in the background PDFs, the signal yields, and the modeling 
of efficiency function (e) are estimated to be 1 — 3%. 

The triple-product for a B meson decay to two vector mesons takes the form 
if ■ {Ti X £2), where if is the momentum of one of the vector mesons. The other two 
vectors and €2 are the polarizations of the two vector mesons. The presence of the 
triple product can be demonstrated by measuring a non-zero value of the asymmetry At = 

r[vt ■ (t^ X t^l) > 0] - r[vt ■{v^xvi)< o]/r[vt ■{v^y<v^) >o] + r[vt ■{v2xvi)< o]. 

In the SM, At^s are very small. Experimentally, the following two T-odd [2, 13 quantities 
are used instead of At, 

^ MA^Al) II ^ lm{A±Al) 

Al + Al + Ay ^~ AI + Al + A] ■ 

The comparison of these triple product asymmetries [A^^ and A\,) with the corresponding 
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TABLE II: The amplitudes obtained from the angular analysis, where the first errors are statistical 
and the second errors are systematic. 



Mode 




(t)K*+ 


Combined 




0.52 ±0.07 ±0.05 


0.49 ±0.13 ±0.05 


0.51 ±0.06 ±0.04 




0.30 ± 0.07 ± 0.03 


0.12;[|:J^ ±0.03 


0.24 ±0.06 ±0.03 


0^9(^11) (rad) 


-2.30 ±0.28 ±0.04 


-2.07 ±0.34 ±0.07 


-2.21 ±0.22 ±0.05 


arg{A_\_) (rad) 


0.64 ± 0.26 ± 0.05 


0.93;[{:|^ ±0.12 


0.72 ±0.21 ±0.06 



30 - 




cos9^. cos9(^ (t)(^ 

FIG. 3: Projections of the transversity angles with results of the fit superimposed. The points with 
error bars show the efficiency corrected data after subtraction of continuum and B K^K~K* 
background. The discontinuity in the fit to cos 9k* is due to the requirement of cos 9k* < 0.8 in 
B+ (j)K*+{K*+ K+7r°) decay. 



quantities for the CP-conjugated decays (A^ and A^.) provide a true T-violating measure- 
ment. 

Additional variables that can be measured through angular analysis are suggested in 
Ref. 01 . By introducing a time-dependent decay rate that can be written as follows 

T{B(B) ^ V1V2) = e-^* (^A. ± Sa. cos(AMt) t Pxa sin(AMt)) gxg. , (4) 

A<o- 

where S, A, and p are expressed as 

Axx = 1{\A,\' + \A,\') , S,, = i(|A,r-|Ar) , (5) 

A^, = -lm{A^A* - , Sx. = -lm{A^A* + A^^A*) , (6) 

A||o = Re{A\\A* + Z,, A^) , Eyo = Re(A|, A* - A,, A^) , (7) 

p±± = lm(^^AlAj^ , pu = -lm(^^A*A}j , (8) 
P±i = Re(^^[AlA, + AtA^]^ , p||o = -Im(^^[Apo±A*A||]^ . (9) 

The subscript i is one of the {0, ||} amplitudes, A^o and A^y are equal to the T- violation 

parameters Aq- = A^ — Aj. and = ^^"^r) respectively. Since the parameters p\a appear 
in terms that are proportional to sin AMt, which vanish in a time integrated measurement, 
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we only report results for the other twelve A and S parameters. The following equations 
should hold in the absence of NP: 

Saa = , S||o = , = . (10) 

Any violation of these relations would be evidence for NP. 

By separating 5-meson and i?-meson samples and a rearranging the fitting parameters 
in the unbinned maximum likelihood fit, we obtain the helicity amplitudes for the B and 
B samples, the triple-product correlations, and other NP-sensitive observables, which are 
given in Table UTTl and HV1 



TABLE III: The measured helicity amplitudes and triple-product correlations in the i?-meson and 
S-meson samples. 



Mode 


i?-meson combined {(f)K*^ and (pK*^ 


) i3-meson combined {4>K*^ and (f>K* ) 


\Ao\' 


0.41 ±0.10 ±0.03 


0.59 ±0.08 ±0.06 


\A±\' 


0.24 ±0.10 ±0.02 


0.26 ±0.08 ±0.04 


arg{A\\) (rad) 


-2.29 ±0.35 ±0.13 


-2.05 ±0.31 ±0.04 


arg{Aj_) (rad) 


0.74 ±0.31 ±0.10 


0.81 ±0.31 ±0.06 




0.21 ±0.08 ±0.03 


0.28 ±0.09 ±0.05 




0.04 ± 0.08 ± 0.02 


0.06 ± 0.06 ± 0.02 



TABLE IV: The observables sensitive to NP extracted from angular analysis. 



A Observables 


S Observables 


Aoo 


0.50 ± 0.06 ± 0.04 


5^00 


= -0.09 ±0.06 ±0.02 


A|||| 


0.25 ±0.06 ±0.02 




= 0.10 ±0.06 ±0.02 


^±± = 


0.25 ± 0.06 ± 0.03 




= -0.01 ±0.06 ±0.02 


A±o (= A^) = 


0.07 ±0.11 ±0.04 




= -0.49 ±0.12 ±0.07 


A±|| (= aI) = 


0.02 ±0.10 ±0.03 




= -0.09 ±0.10 ±0.02 


A||o 


-0.39 ±0.13 ±0.06 


^||0 


= -0.11 ±0.13 ±0.04 



In summary, the decay amplitudes for B — > (j)K* are determined by fitting the angular 
distributions in the transversity basis. The longitudinal polarization component (IAqP) is 
in agreement with Ref. but differs from the naive SM prediction. The measured value of 
shows that both CP-odd (|A_lP) and CP-even (l^oP + l^||P) components are present 
in (j)K* decays with a ratio of 1:3. The phase of A_i_ and Ay differ from zero or — vr by 3.6 and 
3.7 standard deviations respectively. Thus, our data indicate the presence of final state 
interactions. The triple products Aj^ and differ from zero by 2.4(T and 2.8a, respectively 
[l7l |. This may also be an indication of the presence of final state interactions. However, 
the differences between At's and A^'s (= .4^" = A^o,||) that are sensitive to T-violating 
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asymmetry are consistent with zero as well as |5|. The equations, Saa = 0, Sj|o = 0, and 
A_Li = 0, should hold in the absence of NP. Our data does not show any significant violation 
of those relations. 
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